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Abstract 

An exact cubic open string field theory rolling tachyon solution was recently 
found by Kiermaier et. al. and Schnabl. This oscillatory solution has been 
argued to be related by a field redefinition to the simple exponential rolling 
tachyon deformation of boundary conformal theory. In the latter approach, 
the disk partition function takes a simple form. Out of curiosity, we compute 
the disk partition function for an oscillatory tachyon profile, and find that the 
result is nevertheless almost the same. 
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1 Introduction 



Recently there has been remarkable new analytic progress in the study of cubic open 
string field theory (OSFT) [T]. In particular, an exact rolling tachyon solution was 
found [2], related to tachyon matter and decay of an unstable D-brane. The profile 
of the tachyon component of the full string field obtained by [2J from Witten's cubic 
OSFT is 

oo 

T A (x°) = Xe^ x ° + J2(-±) n+1 ^ n f3ne^ nx ° , (1) 

n=2 

where (3 n are positive coefficients^ with a known integral representation. The authors 
of [2] started from the exactly marginal operator 

V = e^7 X , (2) 

constructed the full OSFT solution recursively, adopting the gauge choice of [3], and 
obtained (CQ). Generalizations to superstrings have been reported in pE], and related 
work is also [5]. 

The solution (pQ) has an oscillatory structure, as was suggested to be characteristic 
for the OSFT rolling tachyon by the previous investigations [6J [T] . On the other hand, 
in the boundary conformal field theory (BCFT) description of the same process^, the 
tachyon field rolls monotonously, represented by the simple exponential ([2]). The 
apparent contradiction was addressed in The OSFT string field solution contains 
an infinite tower of other (massive) fields which are sourced by the rolling tachyon 
component. One can perform a field redefinition to boundary string field theory 
(BSFT) [9j§ variables, in such a way that all other fields except the tachyon are 
zero [11]. In the BSFT field coordinatization the tachyon can then turn out to be 
the simple exponential (T2]), while it was oscillatory in the OSFT frame [6J. Thus 
the marginal OSFT solution (where the tachyon component is off-shell) maps to a 
manifestly on-shell form. Further, it maps to the exactly marginal operator which 
gives a BCFT deformation. For the new full OSFT solution of [2] this was shown in 
[T2] . Since the new rolling tachyon solution relates to the known BCFT deformation, 
in particular the time evolution of pressure of the associated tachyon matter has 
already been calculated in [13J, it corresponds to the disk partition function of the 
BCFT with XV ©, 

p(x°) = Z disk (x°) = - - . (3) 

In this note, we are reporting a curious observation. Suppose we were to consider 
BSFT with an oscillatory off-shell tachyon profile of the form (flj). Consider the 



We follow the convention where the true minimum of the tachyon effective potential is at some 
T > while keeping A > 0. We work in units where a 1 = 1. 

2 For another reference on the relation between SFT solutions and deformations of BCFT, see [5]. 
3 A pedagogical discussion of BSFT is also [TO] . 
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worldsheet CFT and turn on the boundary the tachyon field (pQ), 

S = S + I dt T x (X°(t)) , (4) 

JdT, 

it is off-shell and breaks the conformal invariance on the boundary. Suppose we 
attempt to do a straightforward calculation of the disk partition function, leaving the 
zero mode x° unintegrated. Given the oscillatory behaviour of ([1]), we would probably 
expect the resulting disk partition function to be quite unwieldy and very different 
from ([3]). 

However, when we perform the string worldsheet theory analysis (along the lines 
of [13]), surprisingly we find that the result is almost the same as ©, with maximum 
1% relative deviation. The deviation only appears at times close to the value x° ~ 
— ln27rA. Apart from the deviation, there is no oscillatory behaviour - at late times 
the disk partition functions become identical. We do not quite know how to interpret 
this curious observation. Apparently the field redefinitions involved in mapping from 
the oscillatory tachyon profile to the monotonously rolling one are not always so 
significant from the point of view of interesting observables. Further, while in our 
calculation the tachyon is of the form ([1]), the actual values of the coefficients j3 n do 
not matter much - in particular they (and the tachyon field) need not be the same 
as in [2]. Interpretational issues aside, we believe that the calculational tricks which 
we have used will be useful for other investigations and thus interesting in their own 
right. 



2 The disk partition function 

In the first quantized string worldsheet approach, we turn on the tachyon background 
(jl]). The disk partition function is (separating out the zero mode X° = x° + X'° and 
leaving it unintegrated) 



Z disk (x°) = j VX'°VXe- So exp (- j dt T x (x° + X'°(t)) 



(5) 



Note that, in the limit (3 n >i 0, we expect to produce the familiar results for half 
S-brane [T3j . 

By expanding in the boundary perturbation in ([5]) as a power series, and carefully 
following the calculational steps outlined in [14] . the disk partition function is 



^ // i\n\n/3 c nx l '\N rl r , , 

z^(x") = n e (( — /; ' / <--«n e " x ' 0< ' i >> 

n=l N n =Q n ' •* n,i 

= e fn (( " i y )-^.^.-). <6) 

{7Vi,AT 2 ,...}=0 \n=l n ' / 
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with 

and Pi = 1, and where 
I(N U N 2 ,...) 



z n = 2n\ n f3 n e nx > 



(7) 



2tt 



' °° N » jAn) 

nn 

jn=l i=l 
oo N n N m 

n nn 

l<n<m i=l j=l 



n n 

,n=l l<i<j<N n 



\U {n) U {n) \2n 2 

e 1 — e j 



denotes an infinite product of coupled integrals. 

The above formulas are just formal expressions, before good domains of conver- 
gence are found. It is difficult to analyze the problem fully - so we will first study a 
simpler toy model. 



3 A warm-up calculation: the Dyson series 

We have two tasks at hand: (i) to try to calculate the integrals (JSJ) and (ii) to try to 
control the series ([6]). These tasks appear to be rather challenging, so we will first 
consider a toy model calculation. It is reminiscent of the actual one but allows us to 
carry out both tasks. 

We consider a series expansion, which we will call the "Dyson series" from now on. 
It is inspired by the integration formula to compute the canonical partition function 
of a Dyson gas [15], 



N Hi 

i=l 



n 



e iU _ e itj 1/3 



[T(l + f )] 



for which various proofs have been presented in the literature (see [H]). The integral 
§8§ resembles an infinite product of decoupled Dyson gas integrals (El), except for the 
last cross coupling term in the square brackets in the integrand of (jSJ). Let us first 
truncate the infinite product and keep just n max first terms, with integer n max 3> 1. 
(In the end we will consider the limit n max — > oo.) Then, consider the cross coupling 
term in the integrand of ([8]), which renders the integral difficult to evaluate. Let us 
rewrite it as 

r> M IV n N N / -i( m ) \ nm / -An) \ nm 

n nni^"-^ > r m - n nn >-^ i'- e - 

l<n<m i=l j=l l<n<m i=l j=l \ e ! / \ e J 

(10) 

Now it turns out that the integral simplifies drastically if we replace the exponent 
nm in the first term on r.h.s. by n 2 , and the second exponent nm by m 2 . This step is 
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clearly ad hoc. However, it is a useful trick to try, since it simplifies the calculations 
enough to give a tractable toy model calculation to practice with and to gain insight 
for the actual disk partition function calculation. So we consider a version of the 
series ([6]) where we replace the original integrals (jHJ) by 

I{N l ,N 2 ,N^...-n max ) 



where the last line is the exact analytical result for the integral [HH [T7] . Since the 
integrals (fTTT) are a variation of the Dyson gas integral formula (Q, we call the new 
series "Dyson series". In Appendix A we compare the original integrals / with the 
approximate ones /, for some cases where it is possible to calculate the original integral 
analytically, to see how much Dyson series toy model deviates from (iTTj) the exact 
formula. 

The virtue of the Dyson series is that we can also solve the task (ii): we can 
actually sum the series in a controlled way. We will first recognize it as an asymptotic 
series, but can rewrite it as an integral formula which we can regulate by a suitable 
deformation of integration contour. We will discuss that next. 



"Umax N n lAn) 
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4 Summing the Dyson series 

Instead of the series ([6]) we consider the Dyson series with coefficients / instead of /. 
We also simplified it further by truncating the infinite product, so that we have 

^DysonlX , ll max ) - I 11 AT | I TJ^max Tp/i , n 2\]N„ " ^> 

\n=l N„=0 n ' J 11 «=1 [ V n 

Even after truncating to a finite product of n max terms, the expression is not well 
behaved since the product is that of possibly divergent infinite series. In order to gain 
better control, we rewrite ( TT2l) as an integral representation, 

(Umax OO //_-|\ n \N n \ fOO 
n=l N n =0 n ' V ; ' / J0 

POO 

= du exp 
Jo 



Umax 
71=1 
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Now we have a single integral, and the exponent in the integrand is a finite sum of 
n m „~ terms. Let us take a closer look at it. We denote 



u max 

,2 



'-max /_-, \ n „<■ 

F nma Au) = -u+j: { l Jf . (14) 

n=l ^ ' 

For real u, F nmax (u) is oscillatory with the amplitude of oscillation increasing with u. 
The largest oscillations are due to the terms with n ~ n max . As a consequence, the 
integral ffT3l) does not have the expansion ffl2|) for small z n , and the limit n max — > oo 
does not exist. We will next give a prescription to regulate the integral. 

Let us deform the contour of integration in ( fl3i) away from the positive real axis. 
If the integrand would be analytic, this would have no effect. However, it has an 
essential singularity at infinity. Consequently, the contour deformation will change 
the integral, due to a different approach to the point at infinity. Thus we can regulate 
the integral ( Tl~3l) by finding a suitable contour deformation. However, the integral 
will then also become complex valued. Since the pressure is real valued, we adopt a 
prescription where we define it to be the real part of the integral over the deformed 
contoui0, 

Z Dyson (x ;n max ) = Re J du exp [F Umax (u)] , (15) 

where C runs from to oo such that ~ReF nmax decreases monotonically on it. For the 
choice of C, see Fig. 1 which depicts the eye-appealing structure of the real part of 
F nmax (the plot is shown for the value n max = 11). The regular structure of ReF nmax 
arises from the fact that ReF nmax (u) is dominated by the nth term of the sum at 
\u\ ~ n 2 . Fig. 1 suggests that there is a preferred choice for a path (in the quadrant 
< (f) < 7r/2) from to oo that avoids all the light gray regions and proceeds in the 
direction of darker color (decreasing ReF nmax ). We call such a path C pre f and focus 
on ffl5|) with C = C pre f which stays well defined in the limit n max — > oo. 

As an example, let us consider the leading correction with n max = 2. We take 
C pre f with a constant phase, i.e., u = re' ln ^ with r = . . . oo. Then 

POO 

Z Dy son(x°; n max = 2) = Re / dr exp [zvr/4 - (1 + Zl )re i7T/4 - z 2 r 4 /2A] , (16) 

which is well defined. (Recall that z n = z n (x°) ~ exp(nx ).) Developing the integrand 
at z 2 = we get back the (asymptotic) series 



4 With this prescription, it reproduces the asymptotic series (TT^l) . If one has a strong preference 
to keep the integral real valued, one can alternatively first write it as a sum of two identical terms, 
then deform the contour in two opposite ways as mirror images of each other so that the two terms 
become complex conjugates. 
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Figure 1: ReFn(r 2 e^) for r = . . . 12.5 (horizontal axis) and for = 0...n/2 
(vertical axis) with x° = 0. ReFn is small in the dark regions. 



We want to compare this to the leading term 

1 



■> Dyson 



(*°; i) = 



1 + zx 1 + 27c\e x 
Numerical integration of (fl6|) verifies that the total correction with n 



^Dyson (•£ j 2) 



,, , n \X , 1 



y Dyson 



*°; 1) 



10" 



max — 2 is small: 
(19) 



max,x°£R 



and well described by the first few terms of the asymptotic series. Note then that 
at late times the first subleading term is ~ ^^f 5 ~ e~ 3x ° , which is much smaller 
than the leading term ~ z^ 1 ~ e~ x ° . One can argue that at late times x° — >• oo all 
subleading terms are negligible compared to the leading e~ x ° behavior. Similarly one 
finds that the n mnr = 3 correction is even smaller 



■^Dyson(^°; 3) — Z-[) yson (x ; 2) 



J Dyson 



>°; i) 



10 



-7 



(20) 



- max,x £R 



Refining the approximation to larger values of n max produces even more negligible 
corrections. Thus the total correction to the leading result ( fl8l) is at most ~ 10~ 3 in 
the Dyson series, even when n max — > oo. Thus, the approximate result for the disk 
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partition function decays exponentially at late times, 

^DysonO ) = lim Z Dyson (x° , n max ) ~ e~ x ° ; x° -> oo . (21) 

'Fornax >OG 

We will now return back to our original problem, the disk partition function 
The lesson from the Dyson series toy model is that it is useful to truncate the infinite 
products by introducing a 'cut-off' n max and then try to see how much the time 
dependence is corrected as n max is increased. If the additional corrections are more 
and more sub leading, they can be ignored in the limit n max — ► oo. The full series is in 
fact well approximated by just the leading terms as x° — > oo. The partition function 
([6]) turns out to have a similar behavior. 



5 The original disk partition function at late times 

Consider again the exact series ([6]). In our toy model the relevant late-time corrections 
are produced by the first terms in the asymptotic series f|T2|) . It turns out that the 
first terms of the exact series (J5J) can also be calculated analytically, without using 
any approximation for /. The first correction terms are those, where most of the 
N 2 , N 3 , . . . are zero. We denote the integral coefficients of these by 



I n (N u N n ) = I(N U 0, 0, . . . , 0, N n , 0, 0, . . .) 



(22) 



so, e.g., J 2 (iV, 4) = I(Ni = N, N 2 = 4,0,0,...). It turns out we can evaluate the 
integrals 



In(N, 1) 



d An) N (1) 

al l TT al i TT 

2vr 11 2vr 11 

l—l i<j 



e n i — e i 



n 



\ei — e 1 



2/i 



(23) 



This is a well-known Selberg integral, and has previously been applied in the context 
of rolling tachyons in [14] . The result reads 



N 



n-1 



3=0 



j\ (N + n + j)\ 
(n + j)\ (N + j)\ 



(24) 



In particular we find 
h(N,l) N + 2 (N + 3)! 



12 



N + 4 
4 



+ 



iV + 3 
4 



h(N,l) 



1 
4! 
1 
9! 



(N + 4)! (N + 3)! 
iV! + (N- 1)! 



(25) 



(N + 7)\ 
+ 20^- — -b + 10 



(A^ + 6)! + (JV + 5)! 



(JV-1)! (iV-2)! (iV-3)! (JV-4)! 
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where the first terms of the sums are the same as in the Dyson series toy model. 
Thus, we find the corrections to Z disk (eqn. (jUj)) that are linear in z 2 ^. 



■Zdisk (x 



N=0 

1 



l + z 2 



^3 77i r 



+ 



^2(1 - 



AH AH 
2 3 (1 - 10zi + 20^ - lOz? + zf) 



1 + Z a 1 + 



1 + 21 



1 10 



+ ■■• .(26) 



From ([21]) it follows that all higher order linear corrections (those depending on z n 
with n > 4) have similar structures. 

Note that the size of the corrections is slightly larger as in the Dyson series. In 
the latter, at late times the correction linear in z 2 was ~ z 2 z± 5 ~ e~ 3x ° but now we 
find ~ z 2 z~ 4 ~ e~ 2x ° . The correction linear in z 3 is subleading, we find at late times 



— A — 2r l 

z 2 z ~ e 
z 3 z 1 ~ e 6X . 



Moving to higher order, the coefficients I 2 (N,2) apparently also have a formula 
similar to (12511 . We find 



2\N\ 



35 



(N 



+ 77 ^ + 7 )' + 07 (^ + 6)' + + 5)! 



(27) 



AH ' " (N — 1)! (iV-2)! (iV-3)! 
whence the correction to the disk partition function that is quadratic in z 2 becomes 

z|(35 - 77zi + 27z? - z\) 



(2? 



Interestingly, at late times this is of the same order as the linear correction, namely 
~ z 2 z^ 6 ~ e~ 2x ° . As we will discuss below, at the order z 2 we will similarly find 



~ 4 ^i 



■2n-2 



~ e 



-2x u 



and generalizing to order z 3 we will find ~ z 3 z ± 



n 3n— 3 



~ e 



-3x° 



All these are small corrections compared to the leading ~ e~ x decay. 

The above are still a tiny subset of all possible terms in the series (Q, containing 
all possible combinations of monomials of Zi, z 2 , £3, • • •■ But we can estimate their 
late time behavior too. 

Equations (I2"5]h (|2"7|) show that the integers I n (N, 1) and I 2 (N, 2) can be expressed 
as finite sums over binomial coefficients. Using methods outlined in appendix B, we 
evaluated 



I(Ni,N 2 ,...) 



I(N!,N 2 ,...) 



(29) 



U n NJ 

for almost all fixed values of N n for which / < 10 19 . Using these results we then 
discovered a generalizaton of the formulae (1251) . ( 1271) for more complicated sets of 
N 2 ,N 3 ,.... We find that for any Ni = N with fixed N 2 ,N 3 ,..., N nmax (with N nmax > 
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and = N„ . i = N n +2 = • • • ), the I can be written as a finite sum 

1 (N 4- 9-/1' 

/(Ar 1 = jv,Ar 3 ,jv 3 ,...,AU,) = grE c ' (iv-Qi 



£=0 



where S 1 = ^^=2* n2 ^n- The relevant fact for the moment is that the coefficients Ce 
turn out to be independent^] of N. We will give an explicit formula for £ max below. 
The corresponding correction term to Zdisk then becomes 



DC 



6Z disk = H [(-l) n z n ] Nn J2(- z J Ni ( N ' ^3, • • • , N n . 

n=2 N=0 

'N + S- 



ma :c > 



S 



n=2 e=o N=0 

n K-1)%] W " E^nrfprr- 

n=2 e=o 1 + lj 



Importantly, for £ ma2 . we founqj an explicit formula 

E W n - X ) N n\ ~ U ma X + 1 • (32) 



f^vnax 

v — 

-max 



n=2 



The combination of (I32p and the schematic formula (13"Tj) allows us to estimate the 
leading late time dependence of all the correction terms to Zdi S k(^ )- At late times 
the leading part of the generic monomial correction fl3~Tl) is given by the term with 
the highest exponent of z±, i.e., the £ = £ ma x term. Then, combining the late time 
dependences 

Umax TT-max 

U z n n ~ expKE^yi 

n=2 

Ylmax 

exp[(E n ( n ~ i)^)^ - (n max - l)x°) 



n=2 n=2 

than . i 



*1 



n=2 



ii 



in a :v 



z; (S+1) ~ exp[-(E^n)-x°] , (33) 

n=2 



5 The formula ([50]) has been evaluated and verified explicitly (with explicit coefficients Ci), e.g., 
for (7V 2 , N 3 , N4) — (1,1,0), (2, 1, 0), (0, 2, 0) and (1, 0, 1) in addition to the cases discussed above. 

6 Using (pM|) it is straightforward to determine i m ax for the corrections which are linear in z n 
(with arbitrary n = n max ). The general formula (|32p was found by first making an educated guess 
and then testing it with computer calculations. So far we have explicitly verified it up to n max — 4 
but have not yet been able to construct a general proof. 
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we find that the correction term (!3T!) behaves as 



5Z, 



disk ~ e 



(34) 



at late times x° — > +00 and is thus subleading. Thus the leading correction is at 
most of the order e~ 2x ° . 



6 Summary 



We have calculated the disk partition function with an oscillatory tachyon field profile 
([Q) instead of the exactly marginal deformation ([2]). The largest deviations, that we 
have found, from the disk partition function ([3]) of the latter are surprisingly small, 
given by (|26|) and (|28|) . Including the largest one (linear in z 2 ) the disk partition 
function reads 



1 



1 + 2n\e x ° 

1 _ ft 

=0 + 



^2(1 ~ Zl) 
(l + ^i) 5 
fe 23;0 — e 3x °) 



1 + e 



2tt (1 + e 5 



(35) 



where x° = x°+ln 2ir\. Fig. [2] shows the disk partition function with A = 1 and with a 
large value of ft — 15 for better visualization. All the deviations seem to contribute 
around x° = — In 2n\ and become smaller in size. We find the result surprising: 
the disk partition function is very similar to (j3J) although the tachyon profile ([T]) is 
oscillatory and very different from the monotonously rolling (j2J). In particular, the 
oscillatory behavior is almost washed out. 




Figure 2: The disk partition function (!35j) as a function of time t. Here A = 1 and 
we used a large value ~ 15 for fi 2 - For reference, the dashed line represents (J3]). 
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APPENDIX A: A simple method for evaluating / 

Let us study an integral of the form 



J n 



n| n 



I e iti _ £ itj I 2k ij 



(36) 



where kij are integers. This form is a generalization of (jSJ), where the exponents n 2 
and nm are allowed to take any values. The integral may be expressed as a finite 
sum by doing a Fourier transform. The "propagator" from to tj reads 



S (tj -U) = \l- e^-^l 2 ^ = V ( ) 

By inserting this to (|36|) and by doing the t integrals we have 



ij \jnij(tj-ti) 



(37) 



J n 



n e 



m / i— 1 m 

IIm E n *- E ru 

i=l \i=l j=i+l 



(3f 



where only m—1 of the conditions in the (Kronecker) delta functions are independent. 
They can be used to fix the values of n\2, rixs, ... so that 



J n 



n e (-^ 

l<i<j<m mj=—kij 



2k 



n v /. _ v-'- j -.. v" 



n 



j' 



(39) 

This formula can be used to evaluate / for small n and N n . E.g., ^(2, 2) is found by 
letting m = 4, fci2 = 1, ki 3 = ku = ^23 = ^24 = 2, k 3i = 4. Some values are tabulated 
in Table [TJ Note that 

(Ni + 4N 2 )\ 



h(N u N 2 



<I(N U N 2 ) . 



(40) 



In Appendix B we present a more efficient method of evaluating I. 
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Table 1: Comparison of I 2 and I 2 . 





N 2 




h{N 1 ,N 2 ) 


point of interest 








1 


1 




k 


1 

2 
3 
4 



k 
1 
1 
1 
1 


fc! 

(4fc)! 

41 k 

JL '- R 
212 " 

5! _ 5!3! _ 40 
3 — 2! — ^ U 
5 3 . 3 . 2 2 = 5|! = || = 3qq 

2 3 3 2 -5-7 = 7! = 2520 


k\ 

(4fe)! 
4lfc 

5 

30 
210 
1680 




1 

2 


2 
2 


7 2 2 4 = SS = 784 
5-3 3 - 17-2 2 = 9180 


630 
6300 


1(1,2) = [0 + ^/(1,1) 


3 


2 


2 4 ■ 3 • 2371 = 113808 


69300 




1 


3 


3 ■ 2 5 7 2 11 2 = 569184 


450450 


Al,3) = [( 1 4 3 )+H]/(1,2) 


2 


3 


2 6 3 3 7- 13 • 61 = 9592128 


6306300 




1 


4 


2 10 3 2 5 2 7 2 11 2 = 1366041600 


1072071000 


/(1,4) = [( 17 )+2 2 5]I(1,3) 



APPENDIX B: A formula for / using matrix determinants 

In this appendix the integral 



I(N U N 2 ,N 3 ...) 



nn 



2tt 



n n 

n=l l<i<j<N n 



, e ^ n) _ e it ( ™ ) \2nm 



(41) 



_n=l i=l 

oo N„ N„ 

n nn 

_l<n<m i=l j=l 

is transformed to a finite sum over certain integer valued functions. This sum can 
then be used to evaluate / exactly for a given set of {N n }. 
For n = 1 (i.e., = N 2 = N 3 = ■ ■ ■) flU]) becomes 

^= II l^-e^'l 2 . (42) 

Here the integrand is the absolute value squared of the Vandermonde determinant 

2 



|AC*i, 



^v)P 



l<i<j<N 



W 



■ Z 



Jjv-1 
N 



N 



E(-i) n n^ n(feM ( 43 ) 

n fc=i 

where Zfc = exp(itk) and II denotes permutations of 1,2, ... ,N. It is easy to check 
that if (143]) is expressed as a polynomial of {zk}, the constant term in the polynomial 
is equal to Ijv- 
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The Vandermonde approach can be generalized for n > 1 using confluent Vander- 
monde matrices. This can be done by differentiation. For example, 



n 

l<i<j<N 



Zi Zj | 



d 



dz 



N+l 



A(zi, . . . , z N ,z N+1 )\ i 



zn+i=zn 



(44) 



fc jl-ijv+l^l " " - *N \ l N+l L ) Z N 

W 

where ni N = 2 and all other rii = 1. This is the determinant of a confluent Vander- 
monde matrix. 

Generalizing to higher n and N n (with M = ^ n nN n < oo) the integrand in the 
definition of I becomes 



n 



pairs 



(n) _ (m) 
i 3 



2nm 



= \detA\ 



(45) 



where A is the M x M confluent Vandermonde matrix 



d 



s-l 



» 



3-1 



(46) 



The relation between n, k, s and i is (uniquely) determined by 1 < n, 1 < k < N n , 
1 < s < n and £(n, k) + s — % with £(n, k) = Y^n=i m ^m + (k — l)n. The result 
evaluates to 



n 



pairs 



_(«) _ r ( m ) 



2nm 



E e *H«nn 

W 



n fe=l 



9 



s-l 



S ( a -i)! v^r. 



( n )\ »<(n,fc)+s- 1 



(47) 



E £ ™nn 



n k=l 



n 



s=l 



(^(n,fe)+s -!)••• (^(n,fc)+s - S + 1) / ( n )\^(n,fc)+ a -s 



(s-l)! 



(4"')' 



E £ ™nn 



fc=l 



(n) 



E™=1 ^(n,fc) + s 



where = exp(itj^), |z^| 2 = 1 was used in the last step, and the Vandermonde 
matrices in the last form are obtained after antisymmetrization. Note that the com- 
plicated expression £{n, k) is only needed for the pick up the permutation variable % 
with the correct index. 
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The constant term is 

i{n u n 2 ,...) = £ h-i M e jv .. jM n n — — -2 

xA(i^(„ )fc ) + i, . . . , i£( n ,k)+n)^(je(n,k)+l, ■ ■ ■ ,ji(n,k)+n) 
X<H^(n,fc)+l + 1" ie(n,k)+n, jt(n,k)+l + ' ' ' + 3l{n,k)+n) (48) 

where 5(i,j) = 5{j is the Kronecker 5-symbol. For n — 1 the 5 restrictions give simply 
ifc — 3k- Using these the result "simplifies" to 

X A (S(^/( n) fc)+i), . . . , S(ie>( n ,k)+n)) A (^(^/(n^J+l), • • • , S (je'(n,k)+n)) 

(n n 
'(n,k)+s 
S=l 8=1 

where K = M — N\, the first sum goes over all increasing injections S : {1, . . . , K} — > 
{1, . . . , M} [so that i<j<& S(i) < S(j)), and f{n, k) = £(n, k) - 

Due to symmetry, one can add the restrictions ie'( n ,k)+i < ie'(n,k+i)+i (f° r all n > 1 
and 1 < k < N n ), and u>( n ,k)+s < k'{n,k)+s+i, je>(n,k)+ s < J>(n,fc)+s+i (for all n > 1, 
fc, and 1 < s < n) and multiply by the ratio of numbers of terms whence the result 
becomes 

I(N U N 2 ,...) 



/(iV l5 iV 2 ,...) 



U n NJ 

e ' ^-^-^ n n rfw _ iv ...in 2 (50) 

S,{i},{j} ra>lfc=l LV U ^ H 

X A {S(iti( n> k)+i), • • • , 5(zV(n,fc)+ n )) A (5(#/(ri,A:)+l)> ■ ■ ■ , S(jt'(n,k)+n)) 



X 



s=l s=l 



where the prime indicates the presence of the above restrictions. In particular, 
UNi,N 2 ) = 7(^,^,0,0,...) 

N 2 

= E ' e ii-*K e ii-ijr II (^(^fc-l) - Sfek)) (S(j 2 k-l) ~ S(j 2k )) 
S,{i},{j} k=l 

x5 (S(i 2k ^) + S(i 2k ), S(j 2k -i) + S(j 2k )) (51) 

where K = 2N 2 and £'(2,k) = 2{k — 1) was inserted. We have written computer 
codes which evaluate I using the formulae (!50|) . (|5"T1) for a given (but arbitrary) set 
of {N n }. 
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